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ABSTRACT 

We quantify the effects of Galactic extinction on the derived luminosity and diameter 
functions and on the density field of a redshift sample. Galaxy magnitudes are more affected 
by extinction than are diameters, although the effect on the latter is more variable from galaxy 
to galaxy, making it more difficult to quantify. We develop a maximum-likelihood approach 
to correct the luminosity function, the diameter function and the density field for extinction 
effects. The effects of random and systematic photometric errors are also investigated. The 
derived density field is robust to both random and systematic magnitude errors as long as 
these are uncorrelated with position on the sky, since biases in the derived selection function 
and number counts tend to cancel one another. Extinction-corrected luminosity and diameter 
functions are derived for several subsamples of the Optical Redshift Survey (ORS). Extinction 
corrections for the diameter-limited subsamples are found to be unreliable, possibly due to 
the superposition of random and systematic errors. The ORS subsamples are combined using 
overall density scaling factors from a full-sky redshift survey of IRAS galaxies, allowing the 
reconstruction of the optical galaxy density field over most of the sky to 8000 km s~^. 

1. Introduction 

Santiago et al. (1995; hereafter Paper I) presented the Optical Redshift Survey (ORS). 
The ORS contains ~ 8500 galaxies with almost complete redshift information, selected from 
the UGC (Nilson 1973), ESO (Lauberts 1982; Lauberts & Valentijn 1989), and ESGC (Corwin 
& Skiff 1995) catalogues. The sample covers most of the sky with |6| > 20° in both hemi- 
spheres, and consists of two largely overlapping subsamples, limited in apparent magnitude 
and diameter, respectively. This sample is used to derive the galaxy density field to cz ^ 8000 
km and, in future papers in this series, to study its dependence on galaxy morphology and 
spectral properties. Its dense sampling and large variety in galaxy content make the ORS ideal 
for this purpose. Nevertheless, since there is no single optical galaxy catalogue covering the 
entire sky that probes significantly beyond the Local Supercluster, the ORS sample selection 
is based on three different catalogues and thus is not uniform over the sky. There may also be 
systematic trends in the magnitude and diameter scales internal to the different galaxy cata- 
logues from which ORS was derived. In addition, the ORS extends to low Galactic latitudes, 
where variable Galactic absorption introduces an additional source of systematic variations in 
the sample selection. Previous studies of the large-scale distribution of galaxies have either 
worked in passbands where Galactic extinction was unimportant (e.g., Strauss et al. 1990) 
or have restricted themselves to regions of the sky in which the measured extinction was low 
(e.g., da Costa et al. 1994; cf.. Appendix A of Strauss et al. 1992b). In this paper, our aim is 
to measure the density field of galaxies in the face of extinction and photometric errors; this 
paper describes our techniques for doing so, and our tests of these methods. 
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Recovering the galaxy density field requires determining the correct selection function 
for the sample. The selection function quantifies the loss of galaxies due to the magnitude or 
diameter limit and to any other selection effect that may be present. Thus, non-uniformities 
in survey sampling must be correctly incorporated by the derived selection function. The 
selection function is also intimately related to the galaxy luminosity and diameter functions 
in magnitude or diameter-limited surveys, respectively. The latter are cosmologically relevant 
quantities; their shape and dependence on environment density or on galaxy intrinsic properties 
have the potential to constrain different scenarios of structure formation (Blumenthal et al. 
1984; White & Frenk 1991; Cole et al. 1994). A precise determination of the luminosity function 
and diameter function at ^ = is necessary when studying the evolution of these functions 
or when comparing observed galaxy number counts with model predictions (Davis, Geller & 
Huchra 1978, Dekel & Rees 1987, Koo & Kron 1992). Previous determinations of the luminosity 
function of optical galaxies have been carried out by Kirshner, Oemler, & Schechter (1979), 
Sandage, Tammann & Yahil (1979), Davis & Huchra (1982), Kirshner et al. (1983), Efstathiou, 
Ellis, & Peterson (1988), de Lapparent, Geller, & Huchra (1989), Santiago & Strauss (1992), 
Loveday et al. (1992, 1995), and Marzke, Geller, & Huchra (1994), among others. Different 
parametric and non-parametric methods for determining the luminosity function and diameter 
function have been developed and applied; these methods are unbiased by the presence of 
inhomogeneities in the galaxy distribution (Turner 1979, Sandage et al. 1979, NicoU & Segal 
1983, Choloniewski 1986, Efstathiou et al. 1988; see Binggeli, Sandage, & Tammann 1988 
for a review). The diameter function has been investigated by Maia & da Costa (1990), 
Lahav, Rowan-Robinson, & Lynden-Bell (1988) and Hudson & Lynden-Bell (1991). Finally, the 
bivariate distribution of diameters and absolute magnitudes has been analyzed by Choloniewski 
(1985) and Sodre & Lahav (1993). 

In this paper we discuss the effects of Galactic absorption and systematic and random 
magnitude or diameter errors on the derived selection function and density field. In §2 we 
quantify the extinction effects by means of Monte-Carlo simulations. A maximum-likelihood 
method is derived to recover an unbiased estimate of the selection function and of the density 
field. The infiuence of other systematic and random errors is also addressed. In §3, we apply 
the methods developed in §2 in order to recover the ORS luminosity and diameter functions. 
Our main conclusions are presented in §4. Throughout this paper, we work in redshift space; 
thus we are not concerned here about corrections from redshift to real space for galaxies in 
redshift surveys. 

2. Extinction and Systematic and Random Errors 

2.1 The Selection Function and Galactic Extinction 

Any magnitude or diameter-limited sample will become sparser at larger redshifts due to 
the increasing loss of galaxies caused by the adopted apparent magnitude or diameter cut-off. 
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This effect is quantified by the selection function, which expresses the fraction of the total 
population of galaxies that are expected to satisfy the sample's selection criterion. In the case 
of a flux-limited sample, the selection function is given by 

r , , ^{L)dL 

where Lminix) = '^T^f'^ fmin is the minimum luminosity necessary for a galaxy at distance r to 
make into the sample, /mm is the limiting flux^ of the sample in ergs~^ cm~^, and ^{L) is the 
luminosity function in the passband in which the survey is done. In principle, the lower limit 
of the integral in the denominator should be zero, but we set it to a small value Ls = Lmin{fs)i 
as it is impossible in practice to quantify the faintest end of the luminosity function. We set 
Tg = 500 km s~-^ in what follows. In the equation above and throughout this paper we assume 
that ^{L) is universal, showing no dependence on location in the Universe. We will test this 
assumption a posteriori in §3. 

Yahil ct al. (1991; hereafter YSDH) adopt a simple parameterized form for the selection 
function, characterized by three numbers, which experience has shown is sufficiently general 
for these purposes (e.g., Santiago & Strauss 1992): 

0(O=(fe) '"(SSj) ^ r>rs ^ (2) 
I 1 r < rg 

Equation (1) implies that the luminosity function can be related to the selection function by a 
derivative, yielding: 

,,,,.c(",.A^)(±y(.,^y\ (3) 



L ~\~ L J \ J \ L 

where a, j3 and L* = 4:7rr'^fmin are free parameters and C is a normalization constant. Thus 
r* is the distance at which galaxies at the survey fiux limit have a luminosity characteristic of 
that of the knee of the luminosity function. The best-fit values for these parameters are given 
below in §3. This form is a generalization of the Schechter (1976) function with one additional 
free parameter (cf., Yahil 1988). We will use the parameterization of Equation (2) in this 
paper, but Equation (1) will have to be generalized to reflect the effects of extinction. Before 
doing so, let us remind ourselves how one would solve for the selection function for the sample 
in the case of no extinction. Following YSDH, we maximize the likelihood A, the product of 



^ In the case of a survey limited in magnitudes to mum, fmin 10 
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probabilities p that each galaxy have its measured luminosity, given its redshift. For a purely 
flux-limited sample, without the effects of extinction, this is: 



N 



A = Ylp{Li I n) 



i 




^{L')dL 




(f>{ri) dr 



max 



{Li) 



(4) 



where N is the total number of galaxies in the sample and r^ax = f{f / fmin) is the maximum 
distance out to which the galaxy could be placed and remain in the sample. The last step follows 
from Equation (1). 

The derived selection function is used for recovering the density field; the sampling se- 
lection effect is compensated by weighting each galaxy by the inverse of the selection function 
at its position, Wi = l/0(rj), i = 1, . . . ,N. Densities can then be computed by summing each 
galaxy's contribution within a given volume and dividing by that volume. We usually express 
densities in units of the global mean density 



computed over the entire sample volume V out to some limiting redshift (see YSDH, Eqns. 13- 



Now let us consider the effects of extinction. Consider a galaxy whose intrinsic apparent 
photographic B magnitude is given by m. If the galaxy lies in the direction (Z, 6), its observed 
magnitude is given by 



where is the extinction in the given direction and passband. If ruobs is the total 

observed magnitude, = 1; extinction simply reduces the surface brightness at each point 
on the galaxy's two dimensional figure by For isophotal magnitudes, however, not 

only does every point in the galaxy become dimmer, but the radius of the isophote shrinks, 
and thus we expect jrn > 1? the exact value depending on the two-dimensional light profile 
of the object (Cameron 1990). For example, the effects of extinction for a spiral disk with an 
exponential light profile of e-folding length 9^ can be calculated analytically. The observed 
isophotal diameter Bobs is related to the intrinsic isophotal diameter by: 




(5) 



15). 



mobs = m + ^rnAB{l,b), 



(6) 



Oobs^O -0.921 OeAB. 



(7) 
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The quantity 7^ as defined in equation (5) is then 

, 2.5 , l-e-^(l+x) l-e"^ 

lm = l + J^ log i_^-.„..^l^,^^^) ^ l_^-.^l^,) (8) 

where x = 9/9^ and Xobs = ^obs/^e ctre the number of scale lengths subtended by the intrinsic 
and observed diameters, respectively. The second equality holds in the limit of small As- 
Notice that the second term on the right-hand side of Equation (8) depends on the profile 
shape through the number of e-folding lengths within 9 and 9obs- 

We can similarly define a quantity 7^ which gives the fractional decrease in isophotal 
diameter with extinction: 

5 , 6* 2.0 

where again, the second equality holds in the limit of small Ab- 

For a given surface brightness profile, the quantities and 7^ depend only on x and 
Ab- They are related to the fractional changes in flux and diameters by the equations 5f/f = 
1 — iQ-O-^^^B-fm aj^^ §0^Q — I _ iQ-o.2AB'yd ^ respectively. Figure 1 shows 7^ and 7^ as a 
function of 9/9e for Ab = 0.3 mag; the dependence on Ab is weak for small and moderate Ab- 
Panel a assumes an exponential disk and panel b refers to a de Vaucouleurs profile, for which 
expressions analogous to Equations (7)-(9) can be derived. In the latter case, 7^ also depends 
weakly on 9e] we assumed a typical value for early-type galaxies. In both panels, the horizontal 
bars indicate the observed range of 9/9e values for disk and bulge galaxies, assuming limiting 
isophotal levels for detection characteristic of the ESO and Palomar plates. Note that 7^ is 
significantly larger than 7^ for any value of 9/9^; fiuxes are more affected by extinction than 
are diameters. However, 7^ is a stronger function of x than is 7^, since 7^ approaches for 
large extinction, while 7^ approaches 1. Thus the approximation we will use of a single value 
of 7 for an entire sample is likely to be less appropriate for diameter-limited samples than for 
magnitude-limited samples. 

For a fiux-limited sample, the effect of Galactic obscuration is to introduce a non-uniform, 
direction-dependent selection effect: the effective magnitude cut-off limit for a galaxy lying 
in the direction {l,b) is given by mum = mim, observed - lmABil,b), where miim,observed = 
const — 2.5 log is the magnitude corresponding to the fiux limit. In terms of the selection 
function, we now have 

where Lmin{r, ^tuAb) = 'iTrr'^ fmin^O^''^^"''^^ is the minimum luminosity necessary for a galaxy 
in this direction to make it into the sample. A similar expression holds for diameter-limited 
samples. 
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In the absence of absorption = 0), Equation (10) becomes the usual definition 

of the selection function, 0(r) (Equation 1), showing no dependence on direction in the sky. 
We call this quantity the intrinsic selection function </>(r), in order to distinguish it from the 
selection function that incorporates extinction, 006s (r, 7^b). The two functions are related by 
the equation 

4>obs{^nAB) = 4>{rio'-^''^n. (11) 

Notice that it is rather than (f)obs that is directly related to ^{L) via Equation (1). On 
the other hand, it is 0ofes that quantifies the extinction selection effects inherent to the sample. 
We thus need 0(r) to derive the luminosity function, and (f)obs to derive an unbiased density 
field. 

What are the errors resulting from not taking Galactic absorption into account when 
deriving the selection function and density field? We answer this question by means of simu- 
lations of a redshift survey. We use a standard Cold Dark Matter (CDM) N-body simulation 
(0 = 1,/?, = 0.5, A = 0) generated by OL and Karl Fisher, and select a random subsample 
of points around an origin chosen to represent the observer, according to an input selection 
function. Our results on the effect of extinction do not depend on the details of the cosmolog- 
ical model or selection function adopted. Fluxes are assigned to each selected point following 
the luminosity function consistent with the chosen intrinsic selection function. A "Galactic" 
coordinate system is set up as well, and each magnitude is attenuated according to the Burstein- 
Heiles (1982) extinction in that direction. The value of 7 used for each point is drawn from a 
Gaussian distribution with mean unity and dispersion a = 0.2. We thus allow for variations 
in the extinction correction from object to object, and uncertainties in the extinction map. 
The sample was then cut at niobs < 14.5. At this point, we wish to search for biases in our 
techniques for correcting for extinction, as opposed to assessing random errors in the derived 
quantities, and thus we generated a random catalogue with over 10,000 points, roughly twice 
that of the ORS sample itself. 

We first compute the selection function for the sample without attempting to apply any 
correction for extinction effects; we make straight use of the method described in YSDH and 
outlined above. The input (p{r) used in this and all subsequent Monte-Carlo experiments in 
this paper is shown in Figure 2a. The ratio of the derived to the input selection function 
is shown in Figure 2b as a dashed line. The monotonic decrease in the ratio is caused by 
the loss of galaxies from the sample due to Galactic absorption. However, this discrepancy is 
due mainly to the mean extinction over the sample; it is not sensitive to the variation of the 
extinction from point to point. The dotted line in Figure 2b is the quantity 0o6s/0(r), where 
0o6s was calculated from Equation (11), using the input selection function and the mean Ab 
taken over all the points in the simulation. This scaled selection function is very similar to the 
one derived directly from the sample. 
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It is the intrinsic selection function 0(r) that we are interested in recovering. We can 
generalize the maximum- likelihood expressions of Equation (4) to take extinction into account: 



N 



i 





0(lOO-27i^^rj) dr 



1 d(f){r) 



(12) 



where N is the total number of galaxies in the sample. 

Equation (12) takes extinction directly into account when fitting for the parameters of 
0(r). However, we do not have a priori knowledge of the values of 7^, z = 1, . . . ,N. This would 
require knowing the surface brightness profiles of all galaxies in the sample. We thus make 
the approximation of a single value of 7 for all galaxies, as if they all had the same surface 
brightness profile. 7 is then treated as an additional parameter to be fitted simultaneously to 
the luminosity (or selection) function. In Figure 2c wc show the result of maximizing this new 
likelihood function for the same Monte-Carlo simulation as in Figure 2b. The derived selection 
function is now in very good agreement with the input one. Given the intrinsic selection 
function, 4>obs{f-, I^b) for each galaxy follows directly from Equation (11). Thus this approach 
allows us to successfully recover the correct luminosity function (free of extinction biases) for a 
redshift sample, even in the presence of moderate variations in the true value of 7 from galaxy 
to galaxy. 

We next tested our ability to measure 7 itself. We carried out a series of simulations in 
which the true mean value of 7 was varied; in each case, the value of 7 for each galaxy followed 
from a Gaussian distribution with a dispersion of 0.2. Figure 3a shows the derived 7 values as 
a function of the mean input 7. The solid points were derived from simulations of size similar 
to that used in Figure 2, containing about 10,000 points. The open symbols were obtained 
from simulations whose sizes are comparable to the ORS subsamples used in §3 and that had 
random magnitude errors (typically of 0.5 mag) added to the extinction effect. Formal error 
bars from the maximum-likelihood fits are plotted; the smaller simulations have larger error 
bars, as expected. The agreement between derived and input 7 values is still good in all cases, 
although there is a small trend towards underestimating large input 7 for the small sample. 
We investigated the origin of this trend by plotting the derived selection functions (normalized 
to the input one) for each ORS-sized simulation in Figure 3b. Contrary to the error free case, 
(j){r) is now shallower than the input function used (ratio greater than unity). This is caused 
by Malmquist bias, as will be explained in §2.3. Notice, however, that the amplitude of the 



8 



bias in 0(r) is insensitive to 7, indicating that the trend in 7 seen in Figure 3a does not lead to 
discrepant best- fit solutions for 0(r). Instead, the systematics in 7 just compensate for similar 
biases in the 0(r) parameters; indeed, the derived density field (using the methods discussed 
in connection with Figure 5a below) shows no systematic bias that grows with 7. We thus 
conclude that our revised likelihood approach recovers the right (f){r) (apart from Malmquist 
bias) even when 7 is significantly different from unity and in the presence of random scatter 
in the magnitudes. 

We now consider the effect Galactic absorption has on the derived density field. We again 
use the Monte-Carlo simulation with input mean 7 = 1. Galaxies are weighted by the inverse 
of the derived selection function, (/>obs('': 7^b): and smoothed with a tophat of radius given 
by the mean interparticle separation at each point. Densities p(r) are measured on a grid of 
points in space, interspaced by the mean interparticle distance at their location. We normalize 
the densities by ni as given by Equation (5), the global mean computed over the whole sample 
volume to 8000 km s~-^, and define deviations from the mean as 5{r) = p(r)/(p) — 1. 

Figure 4 shows how these grid point normalized densities, 5(r), are affected by Galactic 
absorption. In panel a we plot the errors in the density, A5(r) = 6{r) — 5obs{^), as a function 
of Ab- The quantity 5o6s(r) was computed by weighting each galaxy by the inverse of the 
selection function which completely ignores extinction (Figure 2b). We show the mean value 
of AS for different bins in Ab (solid line). There is clearly a trend in the sense that densities 
are underestimated (overestimated) in regions of high (low) absorption. This is a direct con- 
sequence of the insensitivity of the selection function to variations in extinction over the sky. 
In panel b, AS is plotted as a function of redshift. The mean over all points in each distance 
bin is again plotted as a solid line. Also plotted are the mean curves for the points in each 
bin with the 25% highest and lowest Ab values. Notice that AS for the upper quartile tends 
to increase with distance, whereas the lower quartile shows a systematic decrease: the bias in 
the density field caused by extinction increases with distance. This is consistent with the in- 
creasingly larger disagreement between the derived and input selection functions in Figure 2b. 
At V ~ 4000 km s-\ the error in density can be as large as AS ~ 1.0, although the typical 
values are around 0.2. The results shown are for (7) = 1; these effects become larger for larger 
values of 7. On the other hand, there is no systematic bias in the mean density at a given 
distance; at each distance, there will be directions with both high and low Ab values, whose 
densities are underestimated and overestimated, respectively. We thus conclude that ignoring 
extinction tends to cause biases only in local density estimates, while densities averaged over 
the solid angle of a survey are not systematically affected. 

Finally, in panel c, we show how the difference between intrinsic and derived densities 
behaves as a function of the intrinsic density itself. Again, no systematic trend is seen when 
all points are considered. However, for the high and low Ab quartiles, we see errors of about 
20-25% the mean density. 



9 



Figure 5 shows the results of weighting each galaxy by the inverse of ^obsl'^jT^e) as 
derived by simultaneously fitting 0(r) and 7 and then applying Equation (11) to each object. 
No systematic error in the densities is now visible either as a function oi Ab, distance or 
intrinsic density. In particular, the mean curves for the highest and lowest extinction quartiles 
(panels b and c) are both very close to zero. In addition, the scatter in the diagrams has 
decreased substantially. The remaining scatter is almost entirely due to shot-noise; rederiving 
the density field for the Monte-Carlo sample used in Figure 5, but with exact input values of 
the 7i and of the selection function made almost no difference in the scatter. 

In summary, we have found that the bias in the density field caused by extinction is a 
direct consequence of that in the selection function. In particular, by trying to fit a selection 
function which depends only on distance to a sample affected by Galactic absorption, we do not 
incorporate the variability of the extinction as a function of direction. Local density estimates 
are then systematically in error but densities averaged over solid angle at a given distance are 
not affected. By taking extinction explicitly into account when deriving the selection function 
we recover both the correct galaxy luminosity function and density field, even in the presence 
of scatter in the value of 7 (as long as that scatter is uncorrelated with position on the sky). 
Most of the remaining scatter present in the estimated densities is the result of shot-noise. 

2.2 Other Systematic Effects 

As far as their influence on the derived density field is concerned, systematic errors fall 
into two categories: position-dependent and position-independent errors. In the former case, 
without a priori information on the form of the error (such as we have in the form of the 
Burstein-Heiles maps in the case of extinction), one cannot solve or test for it. In the latter 
case, however, a fit to a radial selection function will absorb the systematic effect^, leaving the 
derived density unbiased. 

Consider, for instance, the case in which the observed magnitudes are incorrect by both an 
offset and a scale error (such as might be caused by uncorrected non-linearities in photographic 
plate material): 

niobs = am -f- 6, (13) 

where a and b are constants and m represents the true magnitude. The selection function 
that accommodates this kind of magnitude error is given by Equation (1) with Lmin{r) = 
47rr^/^''j" lO^'"^^/". Notice that both the selection function and the cut-off luminosity L^i^ are 
still functions of distance only. Our maximum likelihood method will thus find the correct 
selection function even in the presence of such systematic effects, as we demonstrate with a 
Monte-Carlo simulation. Observed magnitudes were obtained from Equation (13) with a = 

^ This is true to the extent that the functional form of the selection function to which we 
are fitting is sufficiently general to be able to absorb the systematic effect. 
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1.04, b = —0.3, giving a systematic effect comparable to tliat of extinction. This simulation did 
not include the effects of extinction or random magnitude errors. The trend given by Equation 
(13) is assumed unknown a priori, and as such is not explicitly corrected for when deriving 
4>{r). Figure 6a shows the ratio of derived to input selection functions. The former is steeper 
than the latter due to the loss of galaxies caused by the systematic trend in the magnitudes. 
The corresponding will be then equally biased. The derived density field, however, shows 
only a small bias, both as a function of distance (panel b) and intrinsic density (panel c). The 
residual bias increases if we make the systematic error given by Equation (13) stronger, since 
the adopted parameterization for (p{r) becomes increasingly inadequate. 

On the other hand, consider an error in which the magnitudes in the Northern Hemisphere 
are systematically brighter than are those in the South, by an amount Aq. If we do not know the 
value of Aq but are aware of a systematic offset between the two hemispheres, we can account 
for the resulting variations in sampling by computing the selection function separately for the 
North and South. By restricting the maximum-likelihood fit to galaxies which are subject to 
the same value of A, we are in practice determining 0o6s(i") = 0(rlO°-^'^). The density field 
calculated by combining the two hemispheres, each with their own selection function, will not 
be biased by sample inhomogeneity. However, if the mean density within each hemisphere is 
not equal to the global mean density the normalization applied to the derived densities will be 
wrong. This means that a sample with systematic errors between different hemispheres cannot 
be used to calculate such quantities as the dipole moment of the galaxy distribution (e.g., 
Strauss et al. 1992c). In order to properly normalize the global density field in such a case, one 
needs external information on the density field, such as supplied by the IRAS redshift survey 
(Strauss et al. 1992a; Fisher et al. 1995). This is the approach we use in §3 below. 

2.3 Random Errors 

Random errors in the magnitudes have two effects on the derived density field in a 
magnitude- limited sample: first, galaxies scatter across the magnitude limit mum- Because 
galaxies with true magnitudes fainter than mum outnumber those brighter than mum , the net 
effect is to augment the numbers of galaxies over the case of no errors, an effect which increases 
with distance. Second, the magnitude errors convolve directly with the luminosity function to 
broaden it, thus making the selection function shallower. These two effects work in opposite 
senses, and thus we might hope that they largely cancel one another in the density field. Let 
the distribution of errors e in the fiuxes be g{e,f), where the distribution can depend on / 
itself. Suppose the density fiuctuation at a point r is S{r); the number of galaxies expected 
within a neighborhood dV of r is^: 

nOO nOO 

Nobsir) = [1 + <^(r)] nidV / dL' / g{eJ)^A7rr^f]de; (14) 

Jl' . (r) ^-oo 

^ We have not included the effects of extinction, for simplicity. 
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L' = 47rr^(/ + e) is the observed luminosity, and Airr'^f is the true luminosity. The observed 
selection function 0o6s(r) is given by an integral over observed luminosities, down to L^j^(r), 
of the observed luminosity function. The observed luminosity function is the convolution of 
the true luminosity function with the error distribution function, thus 

/•oo /-oo 

ni(t>obs{r)= dV ^(e,/)$[47rr2/]cJe, (15) 

JL' . (r) J-oo 

which is the same double integral as in Equation (14). Thus the observed density field is given 

by 

1 + Usir) ^ = 1 + -^(r). (16) 

ni(pobs{r)dV 

Therefore the density field is unbiased by random errors in the fluxes. The one assumption we 
made was that the error distribution was independent of position in the sky. 

We demonstrate this with Monte-Carlo simulations. We assume a Gaussian error distri- 
bution in magnitudes with zero mean. In Figure 7a, we plot the derived selection functions for 
two different choices of Gaussian width a. We also plot the result for the case of a depending 
on magnitude: a = 0.2 for m < 13 and increasing linearly with m to u = 0.5 for m > 14.5. 
The selection function becomes increasingly flatter (ratio increasing) with increasing cr, due 
to Malmquist bias. In panels b and c of Figure 7, respectively, we show the mean difference 
between intrinsic and derived densities for each CclSGj clS ct function of density and distance. 
No systematic differences are found as a function of distance. There is a small tendency for 
underdense regions to have their densities overestimated, and overdense regions to have their 
densities underestimated due to the effect of random errors, but this is at most a 10% effect 
for the examples shown. This confirms the robustness of the density field to random errors in 
the magnitudes. Finally, panel d shows the rms error in the density field as a function of the 
density; although the bias does not grow with increasing a, the error certainly does. 

3. ORS Results 

In this section we derive luminosity and diameter functions for ORS galaxies. We also 
discuss the amplitude of extinction effects on the data and its dependence on morphological 
types. Finally, we show the density field projected on several different slices of the local 
Universe. 

The three catalogues from which the ORS sample was defined are based on separate 
plate material and were done by different observers. In anticipation of differences in the 
various photometric systems, we derived the selection functions separately for each catalogue, 
for both magnitude and diameter-limited samples. Luminosity and diameter functions then 
follow automatically from the inferred 0(r). As we saw in the previous sections, these are 
probably biased by unknown systematic and random errors internal to each catalogue from 
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which ORS was extracted. Systematic errors within each catalogue's hmits are difficult to 
quantify. There are indications, for instance, that the magnitudes from Volume I of the CGCG 
(Zwicky et al. 1961-1968) are systematically fainter by A = 0.3 mag than those from the 
other volumes (Kron & Shane 1976, Paturel 1977). Instead of computing a separate selection 
function for this volume, we explicitly correct its magnitudes just as we do for extinction. As 
for random errors, they are estimated to be of approximately Am ~ 0.09 and Am ~ 0.3 for 
ESO and Zwicky magnitudes, respectively (Lauberts & Valentijn 1989, Kron & Shane 1976, 
Huchra 1976). Errors in the diameters are typically 6logD ~ 0.09 (Paturel et al. 1987, 1991; 
de Vaucouleurs et al. 1991). According to the results shown in Figure 7, the associated biases 
in the density field are typically not larger than 10% of the mean density. The quoted errors, 
however, may depend on such parameters as inclination, morphological type or magnitudes 
and diameters themselves. Since our main interest is in the density field, which we found to be 
less vulnerable to random and systematic errors than is the luminosity function, we opted not 
to try to apply further corrections to the ORS selection functions besides those for extinction 
and for volume I of CGCG. 

In deriving the selection functions, we correct diameters and magnitudes for discreteness 
efi^ects using the method described by Strauss, Yahil, & Davis (1991). For the ESO-m sample, 
we take the 1' diameter cut-off inherent to this sample into account when deriving the selection 
function. This is done by substituting m.m[rmax{Li),rmax{Di)] for rmax{Li) in Equation (12) 
for the likelihood function. In all fits to (f){r), we assume that distances are proportional to 
redshifts in the Local Group rest frame, following Yahil, Tammann, & Sandage (1977) to 
convert from heliocentric velocities. Only galaxies with r < 8000 km s~^ are included in each 
fit, as the samples become unacceptably sparse at higher redshifts (cf.. Figure 4 of Paper 1). 
Galaxies around rich clusters are placed at each cluster's center, following Table 2 of YSDH. 

Each galaxy is weighted by the inverse of the product of the selection function and the 
mean density in each subsample, rii. When combining samples, we multiply by a further weight 
given by the inverse of the relative density of IRAS galaxies from the 1.2 Jy redshift survey 
(Fisher et al. 1995) within the solid angle of each sample and within a redshift of 8000 km s~-^. 
These relative densities are 1.025 for the UGC sample, 1.233 for the ESGC sample, 0.760 
for the ESO sample limited in diameter, and 0.785 for the ESO sample limited in magnitude 
(which covers a smaller solid angle; see Paper 1). 

The derived 0(r) parameters for the various samples are given in Table 1. They are 
labeled by the parent catalogue from which they were drawn, with a suffix -d or -m for diameter 
and magnitude-limited samples, respectively. We also list the number of galaxies used in the 
maximum-likelihood fit, the derived extinction correction parameter 7 and the global mean 
density for each sample. Formal fit errors for the (p{r) parameters typically vary from 10% to 
30%. Errors for 7 are larger, especially for the diameter-limited subsamples, where they can 
be larger than 100%! This is at least partially due to the fact that the approximation that 



13 



7d is the same for all galaxies in the sample is not very good (cf., Figure 1). Another factor 
may be systematic and random errors in the diameters, which may be of larger amplitude 
than the extinction effect. Of course, there is also substantial covariance between the values 
of the selection function parameters and 7. The ESGC-d and UGC-d subsamples have 7 < 
(extinction increases the apparent size of the galaxy) , which is unphysical. 

The errors in 7^, for the ESO and UGC samples are 0.28 and 0.35, respectively, consistent 
with the error bars in the upper panel of Figure 3. The 7 values for ESO-m and UGC-m are 
also consistent with unity, as expected. But how much improvement does the incorporation 
of extinction into the fit to (p{r) cause? The quantity A log A, the logarithm of the ratio of 
the likelihood for the best value of 7 to that with 7 set to zero, is shown in the last column 
of Table 1. The extinction corrections give a large improvement to the likelihood for ESO-m. 
For UGC-m, the extinction corrections may be masked by the larger errors in the Zwicky 
magnitudes. For the diameter-limited samples, the improvement is again small. 

In Figure 8 we show the quantity ni0(r), the number density of galaxies in each subsample 
as a function of distance in the absence of inhomogeneities, for the five ORS subsamples listed 
in the first half of Table 1. Panel a refers to the diameter- limited subsamples, while panel b 
shows the magnitude-limited subsamples. The effect of extinction is quantified by showing, 
in the lower two panels, the ratio of ni0 calculated without and with extinction taken into 
account (the parameters in the former case are not given in the tables). For the magnitude- 
limited subsamples, the (f)obs are steeper than the extinction corrected 0(r), in accordance with 
Figure 2b. For the diameter-limited subsamples, however, only ESO-d shows this behavior. 
Given the larger uncertainties associated with 7^, we decided to fix 7^ = 0.6, a value appropriate 
for ESO and UGC galaxies (Fig. 1). The derived 0(r) parameters for these fits are listed in 
the second half of Table 1. 

Luminosity and diameter functions are plotted in Figure 9. The units of the luminosity 
function are galaxies Mpc~^mag~^, while for the diameter function they are galaxies Mpc~^ (5 
log D/kpc)"-*^. Panel a corresponds to the diameter-limited samples and panel b to the 
magnitude-limited ones. As mentioned earlier, these functions are likely to be affected by 
Malmquist bias. Correcting for it would work in the opposite sense as extinction corrections, 
making the luminosity functions steeper, but unfortunately we do not have a sufficiently ac- 
curate model of the magnitude and diameter error distribution to carry out this correction. 
This makes it difficult to compare the results presented here to previously published luminosity 
functions in the literature. 

The ESO, ESGC and UGC diameter functions are systematically displaced from one 
another. This is due to two effects: the fact that the ESO volume is underdense relative to the 
UGC (ESGC) volume by 30% (60%), as given by the IRAS redshift survey, and the fact that 
the ESO diameters are systematically larger (smaller) than the UGC (ESGC) diameters by 
some 25% (10%) (Paper 1). If we correct for these two effects, all three diameter functions show 
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a much improved agreement. However, the best agreement between ESO and UGC is obtained 
if we assume that Deso = 1-20Z)j/gc5 shghtly smaller a correction than in Paper 1. The 
corrected diameter functions are shown in panel c for all three subsamples. Also shown (solid 
line) is the diameter function of Hudson & Lynden-Bell (1991) derived from the CfAl sample 
(Huchra et al. 1983) using UGC diameters (where we put in the same correction factor of 1.2 
used for UGC-d). Notice that the revised scaling factor between ESO and UGC diameters is in 
closer agreement with Lynden-Bell, Lahav, & Burstein (1989). These latter also matched the 
UGC and ESO diameter functions and found a diameter scaling ratio of 1.17 ± 0.07, but they 
did not have as complete redshift information, and they did not correct for the underdensity 
of the ESO volume. As for the ESGC, its corrected diameter function is a bit steeper than the 
other two, which could be caused by residual non-linearities in the diameters in this sample. 
Finally, as the Hudson &; Lynden-Bell sample overlaps substantially with UGC-d, we are not 
surprised to find good agreement. 

The agreement between the ESO and UGC luminosity functions is very good at bright 
magnitudes. UGC is slightly steeper on the most luminous end. However, these plots do not 
correct for the underdensity of ESO relative to the UGC; with this correction made, in panel d, 
one can match the two luminosity functions above the knee in the luminosity functions if one 
assumes that the ESO magnitudes are systematically brighter by 0.2 mag than those of the 
UGC. However, the agreement at the low-luminosity remains poor. Again, this may be a sign 
of residual non-linearities in the magnitudes in one of the two samples, or perhaps differences 
in the corrections needed for Malmquist bias. In panel d, we also show the luminosity function 
derived by Loveday et al. (1992) from a sample with essentially no overlap with ours. The 
agreement with UGC is very good over the entire range of magnitudes, but we again had to 
apply an empirical shift of 0.2 mag to the Loveday et al. data. 

We now discuss the quality of our best fit solutions for the 5 ORS subsamples considered. 
As discussed by YSDH, the fit quality can be assessed by comparing the observed distribution 
of luminosities to that predicted from the derived ^{L); for each galaxy at a given distance r, 
the predicted distribution function of luminosity is given by Equation (12). Moreover, if ^{L) 
is in fact universal (which is an assumption inherent to the maximum-likelihood technique 
used) and the fitted solution is reasonable, the predicted distribution of luminosities should 
match that observed for any subset of the sample without explicit selection on luminosity. 

The histograms in Figure 10 are the observed luminosity and diameter distributions 
for the 5 ORS subsamples considered in this paper. The smooth curves are the expected 
distributions, given by Equation (12) and the observed distribution of distances. 

We quantify the difference between the observed and predicted distributions by computing 
the statistic: 
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where Nghg and Npred ctre the observed and predicted number of galaxies at each bin in absolute 
magnitude (or diameter), respectively, dof is the number of degrees of freedom available in the 

1 /2 

comparison, and agbs = -^o6s Poisson noise in the observed counts. As in YSDH, only- 

bins with more than 5 objects were used in the determination. This quantity is shown in 
each panel of Figure 10, and is repeated in Table 2. The ratios are close to unity, indicating a 
good quality fit. 

In the remaining lines of Table 2 we show the x^/dof for different subsets of the data, 
defined according to distance, density, extinction amplitude and morphological type. The 
values of x^/dof are in general consistent with unity. There are some exceptions, however. 
The large x^/dof for the lowest As range in ESO-m, for instance, could be indicative of 
problems in the BH extinction maps in the Southern hemisphere. However, we would then 
expect a similarly poor fit in ESO-d, which is not the case. The x^/dof in the last distance 
bin of UGC-m could be due to Malmquist bias. 

There are also several large x^/dof values for some ranges in density which might indicate 
a dependence of ^{L) on environment density. In fact, variations in have been found as a 
function of local density, usually quantified by a change in the correlation function of galaxies 
as a function of luminosity (Maurogordato & Lachieze-Rey 1987; Hamilton 1988; Davis et al. 
1988; Bouchet et al. 1993; Park et al. 1994). This issue is, however, still subject to debate 
(Alimi, Valls-Gabaud, & Blanchard 1988, Valls-Gabaud, Alimi, & Blanchard 1989). Finally, 
several morphological subsamples show large values of x^/dof, indicating that the luminosity 
function may be a function of morphology. Again, dependence of $(iv) on morphology has been 
observed by several authors (e.g., Binggeli et al. 1988, Ferguson & Sandage 1991, Santiago 
& Strauss 1992). Notice, however, that no systematic trends in the quality of the fits are 
found as a function of either extinction, density, velocity or Hubble type, for both diameter 
and magnitude-limited subsamples. Thus, we may conclude that the fits obtained to each 
subsample are reasonable; moreover, our data are consistent with our assumption that the 
luminosity and diameter functions are universal. 

Extinction corrections should also depend on morphology via differences in profile shapes 
as shown in §2.2 (cf., Cameron 1990). In order to assess if this causes systematic errors in 
the selection function, we derived separate 0(r) for different subsets defined according to the 
de Vaucouleurs numerical T types. The selection function parameters, number of galaxies, 7 
and global mean densities for several morphological classes are shown in Table 3. We restrict 
ourselves to ESO-m and UGC-m, since we earlier found that extinction corrections for the 
diameter-limited samples are very uncertain. 

The fits are somewhat more uncertain than those of Table 1 due to the smaller number 
of galaxies in each subsample. We again use the x^/dof statistics, listed in the last column 
of Table 3, to quantify the difference between observed and expected distributions. All, with 
the exception of the spirals in UGC-m, are close to unity. Statistical noise is particularly 
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important for the Irregulars/Dwarfs subset (7 < T < 11), where in spite of the good fit, the 
formal errors in the parameters are substantial. Clear differences in the ^{L) parameters are 
seen as a function of type. In particular, 7 is systematically larger for late-type galaxies than 
for early-types. Notice, however, that the errors in 7 are large enough to render most values 
of this parameter consistent with 1. Furthermore, the existence of correlations among the 
parameters used in the maximum-likelihood fit is also clear from Table 3. In particular, (3 and 
are strongly correlated. Figure 11 illustrates this with likelihood contours for the entire 
ESO-m subsample projected on the /3 — plane (panel a) and the 7 — r* plane (panel h). The 
two contours delimit the 68% and 90% confidence level regions in the parameter values. The 
correlation between j3 and r* stands out clearly. 7 is also mildly correlated with r*. 

Finally, we use our best fit solutions to the luminosity function in ESO-m, UGC-m and 
ESGC in order to estimate spatial galaxy densities. In Paper 1 we showed density contours 
on constant redshift shells; in Figure 12 we show planar slices parallel to the principal planes 
in Supergalactic coordinates. Similar plots of the galaxy distribution have been shown by 
Saunders et al. (1991), Strauss et al. (1992a), Fisher et al. (1995) and Strauss & Willick 
(1995) for IRAS-selected galaxies, and Hudson (1993) for optically-selected galaxies. Gaussian 
smoothing was used, with a width given by the mean interparticle spacing of the UGC-m 
sample at each distance. The zone of avoidance (|6| < 20°) is indicated by shading in each 
panel. The calculations were done to a radius of 8000 km s~^, indicated by the dashed circle 
in each panel. 

The Supergalactic Plane (SGZ = 0) is the central panel in the plot. Some of the more 
dramatic features are labeled here: the Virgo cluster (V) is the overdensity near the origin. 
The Hydra-Centaurus (H-C) and Pavo-Indus-Telescopium (P-I-T) superclusters lie just above 
and below the zone of avoidance at SGX = —3000 km s~^. The IRAS galaxy distribution, 
which continues to lower Galactic latitudes, shows these two structures to be contiguous, 
together making up the Great Attractor. At lower SGY is the Sculptor Void (Sc). The Coma 
cluster (Co) hes at SGY = 6500 km s'^, SGX = 0. The Pisces-Perseus (P-P) supercluster 
lies largely in the zone of avoidance, although much of it is apparent at SGX = 4500 km s"-*^, 
SGY = —3000 km s"-*^; it is contiguous with the Cetus supercluster at lower SGY. 

Slices 3000 km above and below the Supergalactic plane are shown in the two panels 
to the right and left of the central panel respectively. The region above the Supergalactic plane 
is dominated by a large void, and the high SGZ extension of the Great Attractor. The most 
prominent structure in the slice at SGZ = —3000 km s~^ is a part of the Great Wall (Geller 
& Huchra 1989) at SGX= - 2000 km s'^, SGY=6000 km s'^. 

Slices at constant SGX are found in the upper three panels. SGX = (center) cuts 
through the Virgo and Coma clusters, while the right and left panels cut through the Pisces- 
Perseus and Hydra-Centaurus superclusters, respectively. Finally, the lower three panels show 
slices at constant SGY. SGY = lies completely in the Galactic plane, and the SGY = 
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—3000 km s ^ panel cuts through the Southern part of the Pisces-Perseus supercluster, and 
the Sculptor Void. 

4. Conclusions 

We have used Monte-Carlo simulations to assess the amplitude of biases caused by differ- 
ent types of systematic and random errors on the derived selection function and density field 
of a galaxy sample. Our results may be summarized as follows: 

1- Galactic extinction causes a non-uniform selection effect across the sky. The selection 
function derived in the usual way (as a function of radial distance only) is incapable of 
accounting for this effect. It at best reflects the mean loss of galaxies averaged over all 
directions. The same conclusion applies to any other systematic position-dependent bias 
present in the magnitudes or diameters used in selecting a sample. 

2- An unbiased selection function can be recovered from a redshift sample in the case of 
extinction, by taking this effect explicitly into account. However, other (unknown) sys- 
tematic as well as random errors on diameters and magnitudes interfere with extinction 
corrections. 

3- The effect of extinction on diameters is both smaller and more variable than that on 
magnitudes. These two factors make it harder to incorporate extinction corrections into 
the selection function derivation for diameter-limited samples than for their magnitude- 
limited counterparts. 

5- The density field is not strongly affected by systematic effects whose amplitudes do not 
vary across the sky, since these can still be properly incorporated into a radial selection 
function. Furthermore, shell densities averaged over solid angle are not biased even in 
the presence of errors as a function of position on the sky, since they depend on radial 
distance only and, as such, involve an averaging of densities taken over the entire sample 
solid angle. 

6- Random errors also corrupt the derived selection function. As in the case of systematic 
errors, however, they significantly affect the density field only if their typical amplitude 
varies across the sky. 

7- Luminosity and diameter functions are systematically biased by errors in the magnitudes 
and diameters, respectively. Unbiased estimates of $(-L) and $(-D) can be recovered only 
if these errors are known and corrected for when deriving the selection function. This 
applies to both uniform and non-uniform errors, either systematic or random. 

Using the maximum-likelihood technique described in §2.2, extinction corrected selection 
and luminosity functions were derived for the ORS magnitude-limited subsamples. These may. 
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however, be biased by systematic and random errors. For the diameter-hmited subsample, 
the effect of extinction was, as expected, harder to quantify, leading to unacceptably large 
uncertainties in the extinction corrections. We thus chose to apply a ftxed "typical" extinction 
correction for the diameter-limited subsamples. 

We derived separate selection functions for the three different catalogues from which 
ORS was drawn. This was done in order to circumvent biases caused by non-uniform sample 
selection among the three parts of the sky. For comparison, we also derived radial selection 
functions that completely ignore extinction. The amplitude of the extinction bias in ESO- 
m and UGC-m is similar to that expected from the Monte-Carlo simulations. The derived 
selection functions were tested by comparing the observed distribution of luminosities in each 
sample to that predicted from the fit itself. This comparison was done for several subsets of 
each ORS subsample. The agreement between observed and expected distributions is good, 
implying consistency with the implicit assumption that ^(L) is universal. The luminosity and 
diameter functions are consistent with recent determinations in the literature. 

Type-dependent selection functions were derived from ESO-m and UGC-m subsamples. 
We find a trend of increasing 7 for later morphological types. However, this result is tentative, 
given the covariance between the selection function parameters and 7. 

Finally, the ORS density field was plotted in planes parallel to the Supergalactic axes. As 
discussed in the previous sections, this density field is fairly unbiased by random and systematic 
errors, apart from position-dependent systematic trends or strong and variable scatter in the 
magnitudes or diameters. In the next paper of this series (Strauss et al. 1995) we will make 
quantitative comparisons of the density field as traced by different morphological subsets of 
the ORS, and well as comparisons of the density field to that of infrared-selected galaxies. 
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Table 1. ORS Selection Function ParamctcrH 



OctllipitJ 




OL 


P 




7 






ESO-m 


2170 


0.40 


6.25 


-1-1-1 

11100 


-1 

1.00 


10.88 


23.2 


EbO-d 


1614 


0.39 


3.19 


6375 


0.60 


6.87 


1.9 


UGC-m 


2965 


0.36 


9.38 


14060 


0.56 


8.08 


0.9 


UGC-d 


1848 


0.33 


2.83 


4075 


-0.32 


8.60 


0.5 


ESGC-d 


1200 


0.47 


5.26 


9500 


-0.34 


18.24 


0.8 


ESO-d 


1614 


0.39 


3.19 


6375 


0.60 


6.88 


1.8 


UGC-d 


1848 


0.33 


2.76 


4115 


0.60 


9.63 


-2.1 


ESGC-d 


1200 


0.50 


5.44 


10565 


0.60 


20.20 


-3.2 



Note: all samples limited to 8000 km s . 

" km 

^ 10~^ Mpc"'^. No global correction of density from IRAS applied. 
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Table 2. x^/dof for the comparison between predicted 



and observed luminosity (or diameter) distributions of ORS subsamples 



Subsample 


ESO-m 


ESO-d 


UGC-m 


UGC-d 


ESGC-d 


All 


33.6/25 


37.0/25 


27.0/27 


27.2/28 


22.5/22 


V < 2000 


24.1/21 


41.1/21 


15.6/22 


24.1/22 


7.6/13 


2000 < V < 4000 


14.5/15 


19.8/14 


17.0/15 


10.6/11 


17.1/12 


4000 < V < 6000 


7.5/9 


5.4/9 


17.9/10 


1.8/10 


11.4/9 


6000 < V < 8000 


9.5/7 


4.6/6 


29.6/8 


4.6/5 


8.4/6 


S < -0.5 


20.9/23 


45.8/24 


21.6/26 


30.5/27 


15.6/20 


- -0.5 < 5 < 0.5 


3.5/6 


2.2/2 


6.8/11 


1.9/6 


2.7/6 


0.5 < 5 < 2.0 


7.5/12 


9.2/8 


19.5/14 


9.8/12 


14.7/7 


S > 2.0 


15.2/15 


11.0/15 


20.8/22 


21.0/21 


20.7/9 


Clusters 


-/- 


-/- 


10.0/15 


20.3/11 


-/- 


Ab S U.Uo 


oi.o/zi 


zz.u/ zz 


QO Q /OK 

dz.o/zo 


01 .1 / Zo 


io.z/ib 


0.05 <Ab< 0.10 


10.8/13 


8.9/12 


11.5/17 


12.8/19 


22.9/17 


0.10 <Ab< 0.25 


14.1/16 


30.3/16 


13.0/20 


17.2/21 


12.4/16 


0.25 <Ab< 1.00 


12.2/15 


13.8/16 


11.9/10 


8.6/12 


7.7/11 


T < 


34.3/22 


13.0/17 


28.0/22 


9.8/19 


30.4/15 


1 < T < 6 


26.5/20 


40.7/19 


26.1/23 


31.9/22 


11.6/15 


7 <T < 11 


25.0/11 


31.0/16 


18.8/19 


26.1/20 


15.3/11 
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Table 3. ORS Selection Function Parameters for Different Morphological Types 



Sample 


^9 


a 


P 




7 




xVdof 








ESO- 


m 








T < -3 


381 


0.38 


7.49 


15540 


0.51 


1.07 


14.2/16 


-2 < T < 


477 


0.35 


8.96 


11880 


1.01 


3.09 


16.8/17 


1 «^ T ^ zL 


OU± 


u.zu 


IzL 71 


1 7990 


1 l^7 


9 9Q 


99 s /9n 




ozz 


n 1 9 

U.LZ 




Duyu 


1 nn 


n 7fi 

U. / D 


1 1 n /9n 

ii.U/ ZAJ 






— u.uo 


9 AO 


loDU 


1 9 


A Qfi 

41. yo 


19 /I 








UGC 


-m 








T < -3 


401 


0.29 


40.07 


35670 


0.42 


0.57 


13.9/18 


-2 < T < 


590 


0.26 


8.81 


11520 


-0.19 


1.71 


19.2/24 


1 < T < 4 


961 


0.00 


2.93 


5190 


1.00 


1.16 


40.9/25 


5 < T < 6 


680 


0.33 


4.85 


8730 


0.85 


1.90 


23.0/25 


7 < T < 11 


215 


0.86 


8.27 


10050 


0.92 


9.45 


15.2/21 



Note: all samples limited to 8000 km s~ . 
°' km 

^ 10~^ hr"^ Mpc~^. No global correction of density from IRAS applied. 
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1- a) 7 values associated with flirxes (7^, dotted line) and diameters (7^, solid line) for objects 
whose light profile is an exponential and which are subject to an extinction of Ab = 0.3 mag. 
7 is shown as a function of the number of e- folding lengths 9/9^. The horizontal bars give the 
typical observed range of e-folding lengths for ORS galaxies. 6J As in a, for a de Vaucouleurs 
profile. In this case, 9e is the de Vaucouleurs radius. 
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2- a) The input selection function 0(r) used for the Monte-Carlo simulations in this paper. 
Note the logarithmic abscissa, b) The effect of Galactic absorption on the selection function. 
The dashed line represents the ratio of the derived to input selection function when extinction 
is completely ignored, in a simulation with (7) = 1 and extinctions drawn from Burstein & 
Heiles (1982). The dotted line is the selection function ratio, where the selection function has 
been scaled using the mean Ab value in Equation (11). Note the linear scale, c) The ratio of 
the selection function derived taking extinction into account, to (f). 
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3- a) Relation between the input and derived values of the extinction correction parameter 
7 from Monte-Carlo simulations. The solid points correspond to samples containing over 
10000 galaxies. The open symbols refer to samples comparable in size to the ORS subsam- 
ples discussed in §3 (about 2000 galaxies). Error bars are the formal uncertainties from the 
maximum- likelihood fit. b) Ratio of derived to input selection function for the ORS-sized 
simulations whose observed and input 7 values are given by the open points in panel a. 
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4- The effect of Galactic absorption on the derived density field. The three panels give the 
difference between the input and derived densities calculated on a grid of points for a Monte- 
Carlo sample with input (7) = 1. The simulated galaxies were weighted by the inverse of the 
derived selection function given in Figure 2b. For clarity, only every second point is plotted, a) 
Density error as function of Ab- The solid line is the mean error in each bin of Ab- b) Density 
error as a function of distance. The solid line is the mean error in each bin of distance. The 
dashed line (dot-dashed line) is the mean error over the points with the 25% lowest (highest) 
Ab values, As in b, now as a function of the input density. The upper and lower envelopes 
in this case are due to the fact that both 6 and Sobs have lower limits of —1. 
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5- As in Figure 4, with extinction properly corrected for in the density calculation. 
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6- a) Ratio of derived to input selection functions for a Monte-Carlo simulation whose magni- 
tudes were systematically biased according to Equation (13) in the text, b) Difference between 
input and derived densities as a function of velocity distance. The solid line is the mean 
difference in distance bins, c) Errors in density plotted as a function of input density. 
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7- a) Ratio of derived to input 0(r) for three Monte-Carlo simulations whose magnitudes were 
subjected to random errors Gaussian distributed with standard deviations 0.2 mag (solid), 
0.5 mag (dashed) and linearly increasing with magnitude (dot-dashed), b) Mean error in the 
density field as a function of velocity for the three Monte-Carlo simulations whose selection 
functions are given in a. c) The mean error for the three simulations as a function of the input 
density, d The rms error for the three simulations as a function of the input density. 
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8- a) Product of derived mean density and selection functions for the diameter-limited sub- 
samples: ESO-d (solid line), UGC-d(dotted line) and ESGC-d (dashed line), b) As in a for 
the magnitude- limited subsamples. c) Ratio of the extinction-uncorrected to the extinction- 
corrected ni(f)(r) values for subsamples shown in panel a. d) Ratio of the extinction-uncorrected 
to the extinction-corrected ni4){r) values for subsamples shown in panel b. 
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9- a) Derived diameter functions for the ESO-d (dashed hne), ESGC-d (dotted hne) and 
UGC-d (dot-dashed hne) subsamples. The units are galaxies Mpc~^ (5 log D/kpc)"-*^. b) Lu- 
minosity functions for the ESO-m (dashed) and UGC-m (dot-dashed) subsamples. The units 
are galaxies Mpc""^ mag~^. c) ESO-d, ESGC-d and UGC-d diameter functions after taking 
into account the differences in mean density and diameter scale between the three samples. 
The solid line is the diameter function of CfAl galaxies by Hudson & Lynden-Bell (1991). 
d) ESO-m and UGC-m luminosity functions after taking into account the differences in mean 
density and magnitude scale between the two samples. The solid line is the luminosity function 
of Loveday et al. (1992), without any correction for Malmquist bias. 
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10- Panel a) The histogram shows the luminosity distribution of galaxies in the ESO-m sub- 
sample. The solid line is the expected distribution from the best- fit solution to ^{L) as given in 
Table 1. The goodness-of-fit parameter x^/dof is also shown. Panel b: as in a for the UGC-m 
subsample. Panel c: as in a for the diameter distribution of the ESO-d subsample. Panel d: 
as in c for the UGC-d subsample. Panel e: as in c for the ESGC-d subsample. 
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11- 68% and 90% confidence level contours in the likelihood function of ESO-m as projected 
on the /? — r* (panel a) and 7 — r* (panel b) planes. 
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